A new upper limit on the total neutrino mass from the 2dF Galaxy Redshift Survey 
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We constrain f v = f2„/Q m , the fractional contribution of neutrinos to the total mass density 
in the Universe, by comparing the power spectrum of fluctuations derived from the 2dF Galaxy 
Redshift Survey with power spectra for models with four components: baryons, cold dark matter, 
massive neutrinos and a cosmological constant. Adding constraints from independent cosmological 
probes we find /„ < 0.13 (at 95% confidence) for a prior of 0.1 < f2 m < 0.5, and assuming the scalar 
spectral index n = 1. This translates to an upper limit on the total neutrino mass m„,tot < 1.8 eV 
for 'concordance' values of f2 m and the Hubble constant. Very similar results are obtained with a 
prior on fi m from Type la supernovae surveys, and with marginalization over n. 



Whether neutrinos are massive or not has been an open 
question for a long time, but the recent data from atmo- 
spheric and solar neutrino experiments (l]-|7| are most 
naturally interpreted in terms of neutrino oscillations, 
which implies that not all neutrinos are massless [p|-^0[. 
However, since the oscillation probability depends on 
the mass-squared differences, and not on the absolute 
masses, the oscillation experiments cannot provide abso- 
lute masses for neutrinos. The mass scale can in principle 
be obtained from e.g. the energy spectrum in the beta- 
decay of 3 H plj , or from neutrinoless double beta decay 
jl2],[l3|]. At the present, cosmological data ]l4|-]l6[ provide 
stronger constraints on the total neutrino mass than par- 
ticle physics experiments. Since neutrinos with masses of 
the order of a few tenths of an electron volt (eV) can have 
a significant effect on the formation of large-scale struc- 
tures in the Universe, observations of the distribution of 
galaxies can provide us with an upper bound on the den- 
sity of massive neutrinos. We will in this paper use data 
from the 2dF Galaxy Redshift Survey (2dFGRS), which 
is the largest existing redshift survey [jl7|,[l8|, to obtain 
an upper bound on the sum of the neutrino masses. 

Massive neutrinos make up part of the dark matter 
in the Universe. In the cosmological model favoured 
by data on large-scale structure and the observed fluc- 



tuations in the Cosmic Microwave Background (CMB) 
|ll9| , [2"ol| , the Universe is flat, and the contributions to 
the mass-energy density in units of the critical density 
are J7a ~ 0.7 from vacuum energy or a 'quintessence' 
field, and f2 m « 0.3 from matter. Baryons make up 
only /b = ^b/^m ~ 0.15 of the matter contribution 
p9| , pfH , most of the remaining being in the form of cold 
dark matter (CDM), the exact nature of which is still 
unknown. 'Cold' in this context means that the par- 
ticles were moving at non-relativistic speeds when they 
dropped out of thermal equilibrium. Particles drop out 
of equilibrium roughly when their interaction rate falls 
below the expansion rate of the Universe. For neutrinos 
with masses in the eV range this happened when they 
were still relativistic, and so they will be a 'hot' compo- 
nent of the dark matter (HDM) . This has implications for 
large-scale structure, since the neutrinos can free-stream 
over large distances and erase small-scale structures (see 
e.g. |2l|] for an overview). As a result, mass fluctuations 
are suppressed at comoving wavenumbers greater than 
k nr = 0.026(m„/l eV) 1 / 2 ^ 2 /iMpc" 1 §§, where m v is 
the neutrino mass of one flavour. The neutrino contribu- 
tion to the total mass-energy density, f2„, in units of the 
critical density needed to close the Universe, is given by 
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where m^tot is the sum of the neutrino mass matrix 
eigenvalues, and the Hubble parameter Ho at the present 
epoch is given in terms of h as Hq = 100 h kms -1 Mpc -1 . 
Eq. ([!]) assumes that all three neutrino flavours drop out 
of equilibrium at the same temperature, which is a rea- 
sonable approximation. The suppression of the matter 
power spectrum P m (k) on small scales is given approxi- 
mately by (see e.g. Q) 
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where /„ = f2„/f2 m . Therefore even a neutrino mass as 
small as 0.1 eV gives a reduction in power of 5-15 %. 

The 2dFGRS has now measured over 220 000 galaxy 
redshifts, with a median redshift of z m « 0.11, and is the 
largest existing galaxy redshift survey jT^. A sample of 
this size allows large-scale structure statistics to be mea- 
sured with very small random errors. An initial estimate 
of the convolved, redshift-space power spectrum of the 
2dFGRS has been determined flf] for a sample of 160 000 
redshifts. On scales 0.02 < k < 0.15 /iMpc -1 , the data 
are robust and the shape of the power spectrum is not 
affected by redshift-space or non-linear effects, though 
the amplitude is increased by redshift-space distortions. 
These data and their associated covariance matrix form 
the basis for our analysis. 

For each model, we calculate its linear-theory mat- 
ter power spectrum, and for the 2dFGRS power spec- 
trum data it is sufficiently accurate to use the fitting 
formulae derived in [E3 . The relation between the mea- 
sured galaxy power spectrum and the calculated matter 
power spectrum is given by the so-called bias parame- 
ter b 2 = P g (k) / P m (k). By definition, b is in principle a 
function of scale, and several biasing scenarios have been 
proposed [E5[. This issue is of some importance for our 
analysis. For example, if the galaxy distribution is more 
biased on small scales than on large scales, a non-zero 
best-fit value for /„ may be obtained. However, on the 
scales we consider there are theoretical reasons to expect 
that b should tend to a constant [M. For the 2dFGRS, 
a recent analysis j27j ] looking for deviations from linear 
biasing found no evidence for it. We will therefore in the 
following assume that the biasing is scale-independent. 
In particular two independent analyses [^,^7]] suggest 
that that the data are consistent with b rj 1 on large 
scales. We choose to avoid the complications in the nor- 
malization of the power spectra caused by redshift-space 
distortions and the actual value of the bias parameter by 
leaving the amplitude of the power spectrum, from here 
on denoted by A, as a free parameter, and then marginal- 
ize over it. 

We shall consider here a model with four components: 
baryons, cold dark matter, massive neutrinos (hot dark 



matter) and a cosmological constant. As an illustration, 
we show in Fig. [I] the power spectra for Q, u = 0, 0.01, 
and 0.05 (all other parameters are fixed at their 'con- 
cordance model' values given in the figure caption), after 
they have been convolved with the 2dFGRS window func- 
tion, and their amplitudes fitted to the 2dFGRS power 
spectrum data. For the 32 data points, the — 0-model 
had x 2 = 32.9, Q v = 0.01 gives x 2 = 33.4, whereas the 
model with f2„ = 0.05 provides a poor fit to the data 
with x 2 = 92.2. 

Clearly, the inference of the neutrino mass depends on 
our assumptions ('priors') on the other parameters. We 
therefore add constraints from other independent cos- 
mological probes. The Hubble parameter has been de- 
termined by the HST Hubble key project to be h = 
0.70 ± 0.07 and Big Bang Nucleosynthesis gives a 
constraint U^h 2 = 0.020 ± 0.002 on the baryon density 
fjOf . For these parameters, we adopt Gaussian priors 
with the standard deviations given above. 

Perhaps the least known prior is the total matter den- 
sity f2 m . The position of the first peak in the CMB power 
spectrum gives a strong indication that the Universe is 
spatially flat, i.e. fl m + fl A = 1 [|l|,|o). The CMB peak 
positions are not sensitive to neutrino masses, because 
the neutrinos were non-relativistic at recombination, and 
hence indistinguishable from cold dark matter. Although 
the shape of the power spectrum is independent of cur- 
vature, the curvature does affect the choice of priors on 
f2 m , and we choose to consider flat models only. When 
the constraint of a flat universe is combined with surveys 
of high redshift Type la supernovae pl]j3"2] | , one finds 
fl m = 0.28 ± 0.14. However, studies of the mass-to- 
light ratio of galaxy clusters find values of £l m as low 
as 0.15 p3[ , whereas cluster abundances give a range of 
values il m ~ 0.3 — 0.9 Another measurement of 

Q m ~ 0.25 which is independent of the power spectrum 
of mass fluctuations and the nature of dark matter has 
been derived from the baryon mass fraction in clusters of 
galaxies, coupled with priors on Q^h 2 and h [ ^Opl]] . We 
will therefore use two different priors on J7 m . The first 
is a Gaussian centered at f2 m = 0.28 with standard devi- 
ation 0.14, motivated by [M. As an alternative, we use 
a uniform ('top hat') prior in the range 0.1 < f2 m < 0.5. 
Given that we use the HST Key Project result J2{| for h, 
Q m < 0.5 is required to be consistent with the age of the 
Universe (42) being greater than 12 Gyr. 

Finally, the CMB data j|||l|4|] are consistent with 
the scalar spectral index of the primordial power spec- 
trum being n = 1, but we also quote results for n = 0.9 
and n — 1.1, and for the case when we marginalize over 
n with a Gaussian prior n = 1.0 ± 0.1, motivated by the 
latest data from VSA and CBI @Q. 

Results will be presented for the case of N v = 3 equal- 
mass neutrinos, but the derived upper bound on the total 
neutrino mass is only marginally different for N v = 1 or 
2. For each set of parameters, we computed the the- 
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oretical matter power spectrum, and obtained the % 2 
for the model given the 2dFGRS power spectrum. We 
then calculated the joint probability distribution func- 
tion for f v and T = tt m h (which represents the shape 
of the CDM power spectrum) by marginalizing over A, h 
and /b weighted by the priors given above. For A we 
used a uniform prior in the interval 0.5 < A < 10, where 
A = 1 corresponds to the normalization of the 'concor- 
dance model', discussed in ]28| ]. Using instead a prior 
uniform in log A, or fixing A at the best- fit value had 
virtually no effect on the results. We evaluated the like- 
lihood on a grid with 0.1 < Q m h < 0.5, < /„ < 0.3, 
< / b < 0.3, 0.4 < h < 0.9, and 0.5 < A < 10. In |l| it 
was found that the 2dFGRS data alone allow a solution 
with a high baryon density /b = 0.4, in addition to a 
low density-low baryon density solution. However, given 
the above priors, the solution with high baryon density 
gets little weight and the fitting-formulae in are suf- 
ficiently accurate for the measured BBN baryon density. 

The results are shown in Fig. || for the cases of no 
prior on Q m (left panel) and with the uniform prior 
0.1 < fl m < 0.5 (right panel). Marginalizing the dis- 
tributions in the right panel of Fig. || over £l m h, we 
get the one-dimensional distribution for f v given by the 
solid line in Fig. ||, and an upper limit f u < 0.13 at 95% 
confidence. For comparison, marginalizing without any 
priors, the limit becomes f v < 0.24. Adding just a prior 
on f2 m , we find /„ < 0.15, so this is clearly the most im- 
portant prior. Marginalizing with just a prior on h or on 
Ob/i 2 , the 95 % confidence limit becomes /„ < 0.20. As 
a further test of the stability of our analysis, we used the 
full set of priors, but only the power spectrum data at 
scales k < 0.1 /iMpc -1 . In this case the limit increases 
to f v < 0.20. 

There is a further degeneracy of f v with the scalar 
spectral index n, since increasing n increases power on 
small scales and leaves more room for suppression by the 
massive neutrinos. Also shown in Fig. |3J are the distri- 
butions for the cases n — 0.9 (dotted line) and n = 1.1 
(dashed line). With n = 1.1, the 95% confidence limit on 
f v increases slightly to 0.16. The results are summarized 
in Fig. H and in Table Q. Also included in the this ta- 
ble are the results obtained using the Type la supernova 
prior, and it is seen that the results for the two different 
choices are almost identical. Running a grid of models 
with n as an added parameter, and marginalizing with a 
prior n = 1.0±0.1, consistent with the CMB data 
we find f v < 0.16 at 95 % confidence. 

To summarize, we have analyzed the shape of the 2dF- 
GRS power spectrum to obtain an upper bound on the 
fractional contribution of massive neutrinos to the total 
mass density, /„, and found an upper limit /„ < 0.13 
at 95% confidence for 0.1 < f2 m < 0.5 and the scalar 
spectral index n — 1. This translates into a constraint 
on the sum of the neutrino mass matrix eigenvalues 
"Vtot < 1-8 eV for Q m h 2 — 0.15. With marginaliza- 



tion over n with a prior n — 1.0 ± 0.1, the limit be- 
comes m„.tot < 2.2 eV. Previous cosmological bounds on 
neutrino masses come from data on galaxy cluster abun- 
dances |l4|,|l| , the Lyman a forest |nj , and compilations 
of data including the CMB, peculiar velocities, and large- 
scale structure |T(| . They give upper bounds on the total 
neutrino mass in the range 3-6 eV. Note that the fluctu- 
ation amplitude derived from the cluster abundances is 
still under some debate p4|-^7|. The most recent limit 
is that of p(J, from a combined analysis of CMB and 
large-scale structure data an upper limit m u ,tot < 4.2 eV 
was found. Our bounds, summarized in Table Q, are 
stronger, largely because of the small statistical errors in 
the 2dFGRS power spectrum, although the priors on the 
other parameters, in particular on f2 m , are also impor- 
tant. Note also that all of these results are stronger than 
current constraints from particle physics. As they stand, 
the controversial results of [|l3| imply a nearly degenerate 
neutrino mass matrix and a bound on neutrino mass ma- 
trix eigenvalue 0.1 < m v < 1 - 20 eV [^6]-^8j, where the 
upper limit is somewhat model-dependent. Our results 
are consistent with this range. If the largest neutrino 
mass is in fact of order a tenth of an eV, it should be 
possible to measure its value using a combination of cos- 
mological data, combining 2dFGRS, SDSS, MAP and/or 
Planck. 
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FIG. 1. Power spectra for Q. v — (solid line), f2„ = 0.01 (dashed line), and Sl v = 0.05 (dot-dashed line) with amplitudes 
fitted to the 2dFGRS power spectrum data (vertical bars) in redshift space. We have fixed Q m = 0.3, S1a = 0.7, h — 0.7, 
flbh 2 = 0.02. The vertical dashed lines limit the range in k used in the fits. 
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FIG. 2. 68 (solid line), 95 (dashed line) and 99% (dotted line) confidence contours in the plane of /„ = Jl„/!l m and F = Q. m h, 
with marginalization over h and Qbh 2 using Gaussian priors, and over A using a uniform prior in 0.5 < A < 10. Fhe left panel 
shows the case of no prior on Q m , and the right panel the case of a uniform 'top hat' prior on fi m in 0.1 < fi m < 0.5. 




0.05 0.1 0.15 0.2 0.25 

f„= aj !l m 

FIG. 3. Probability distributions, normalized so that the area under each curve is equal to one, for /„ with marginalization 
over the other parameters, as explained in the text, for N v = 3 massive neutrinos and n = 0.9 (dotted line), 1.0 (solid line), 
and 1.1 (dashed line). 
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fi m = 0.28 ± 0.14 (SNIa) 
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0.9 


0.12 


1.5 
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TABLE I. Summary of 95% confidence upper bounds on /„ with our two chosen priors on fl m . The conversion of f v to 
nVtot is for h — 0.7 and the central values of f2 m = 0.28 (SNIa case) and Q, m = 0.30 (uniform prior case). 
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